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For a long time a crystal structure of the high- pressure e-phase of solid oxygen was a mystery. 
Basing on the results of recent experiments that have solved this riddle it is shown that the magnetic 
and crystal structure of this phase can be explained by strong exchange interactions of antiferromag- 
netic nature. The singlet state implemented on quaters of O2 molecules has the minimal exchange 
energy if compared to other possible singlet states (dimers, trimers). Magnetoelastic forces that arise 
from the space dependence of the exchange integral give rise to transformation of 4(02) rhombuses 
into the almost regular quadrates. Antiferromagnetic character of exchange interactions stabilizes 
the distortion of crystal lattice in e-phase and impedes such a distortion in the long-range a- and 
5-phases. 
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Solid oxygen is known to occupy a particular place in the family of cryocrystalsi Steady interest to oxygen during 
almost 50 years is due to the magnetic properties of O2 molecule which posses nonzero spin So 2 = 1 m the ground 
electronic state. This ensures magnetism of all the solid O2 phases. 

Solid oxygen has a rather complicated phase diagram, which includes three low temperature phases stable at 
ambient pressure: 7 (T < 43.8 K), (3 (23.8 < T < 43.8 K) and a (T < 23.8 K), and several high pressure phases: 
6, e and £ (see review^ which includes also a comprehensive history of a problem). All these phases except C-O2 are 
insulators. Metallization of solid O2 takes place at 96 GPaM at room temperature. Moreover, C-O2 was also observed 
in a superconducting stated. 

An interesting feature of all (except 7-O2) the phases of solid oxygen is parallel alignment of the molecules which 
is usually explained by strong contribution of exchange interactions into anisotropic (i.e. depending on the mutual 
orientation of molecules) part of intermolecular potential. Presence of stable orientation ordering enables to sim- 
plify substantially many of theoretical models and, in particular, makes it possible to describe the structural phase 
transitions in solid oxygen disregarding orientational dynamics of molecules. 

In the absence of the orientational degrees of freedom, the low-temperature rhombohcdral (space group Rim) (3- 
phase can be thought of as a para-phase for all the magnetic phases. In particular, (3-0% has a planar structure, 
consisting of close packing of parallel oxygen molecules with centers of mass in the apexes of regular triangles, and 
with orientation perpendicular to the basal planes. The temperature dependence of magnetic susceptibility of f3- 
O2 is typical for antiferromagents (AFM)§. Noncollinear 3-sublattice ordering in this phase (Loktev structure) was 
predicted in 7 - and now is generally accepted. Below 23.9 K (or at high pressure) /3-phase becomes unstable and 
transforms into monoclinic (space group C2/m) a-phase. Corresponding a/3-phase transition has magnetoelastic 
nature associated with strong dependence of exchange interaction vs intermolecular distance, as it was shown in-. 
The a-phase possesses the collinear (Neel) magnetic structure with the easy direction parallel to the monoclinic axis 
b of slightly distorted (compared to regular hexagonal) lattice. Long-range AFM ordering, which can be described 
within a simple 2-sublattice model, is stabilized by the "deformation-induced splitting" of intra- and inter-sublattice 
exchange integrals. It should be stressed that all the in-plane exchange constants originate from a single constant 
J(r) that has AFM character (i.e., J(r) > 0), is isotropic and describes intermolecular spin interactions in /3-phase^. 
Mutual shift of the close-packed basal planes that accompanies formation of AFM ordering also has magnetoelastic 
nature and originates from space dependence of inter-plane exchange integrate. 

Hydrostatic pressure up to 4-^6 GPa induces continuous shift of the basal planes, while the magnetic structure of 
a-02 and orientation of molecules remain invariable. At approximately 6.5 GPa the mutual shift of neighboring planes 
attains 1/2 of an in-plane intermolecular distance and solid O2 transforms into orthorhombic (space group Fmmm) 
5-phase. The type of magnetic order in a- and (S-phases is similar (collinear AFM structure) within a6-plane, but 
relative orientation of spins in the neighboring planes (between the first interplane neighbors) is different - parallel 
in a-02 and antiparallel in S-O?^- Due to crucial change of magnetic structure (from 2 to 4-sublattice) a<5-phase 
transition is classified as the 1st order—. Abrupt change of magnetic order at the a<5-transition point is accompanied 
by discontinuous shift of the close-packed planes, which also originates from space dependence of the inter-plane 
exchange integral. 

Further increase of pressure up to 8 GPa produces another phase transformation into e-phase which elusive structure 
has been determined recentl y 12 i 13 . The transition is undoubtedly of the 1st order and is accompanied by a considerable 
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(up to 5.4%) volume reduction. Crystal structure of e-phase is layered, as it is the case for a-, f3-, and <5-phases, and 
has monoclinic (space group C2/m) symmetry. Variation of interplane distance (equal ~ 3.4 Aat the transition point) 
with pressure is very small. So, the volume change is mainly due to the variation of intermolecular distances within 
the basal plane. The peculiar feature of e-phase is association of four O2 molecules into rhomb-shaped (according to 
Refll2l) or square-shaped (according to Refll3l) (02)4 molecular units, which are symmetry equivalent and centered 
on the lattice points at (0,0,0) and (0.5,0.5,0). Common spin-state of the (02)4 = Os cluster is nonmagnetic*^, with 
the total spin So 8 = 0. 

Physical reasons of such an unusual behavior of the magnetic molecular crystal are not yet clearly understood. 
Should e-phase be considered as a chemically new substance, what is the nature of forces that keep (02)4 quadrates in 
neighboring planes locked under high pressure, what is the role of magnetic interactions - all these questions are still 
open. First-principles calculations^ demonstrate the tendency of the O2 molecules for dimerization and formation of 
herringbone-type chains but failed to prove that the (02)4 structure has the lowest energy. 

With the account of these results, in the present paper we make an attempt to elucidate the role of exchange inter- 
actions in formation of nonmagnetic e-phase and show how the pressure-induced variation of the exchange constants 
may produce strong distortion of crystal lattice. 

I. INTUITIVE CONSIDERATIONS 

Analysis of the magnetic and structural properties of /?-, a- and ^-phases of solid oxygen shows^ that exchange 
interactions in this crystal are so strong that they are responsible not only for variation of magnetic order but also 
produce rather noticeable deformations of the crystal lattice. So, it seems reasonable to assume that e-phase makes 
no exclusion and its complicated and surprising structure is mainly due to strong exchange inter-molecule interaction 
that keeps quarters of O2 molecules as the independent chemical units, equalizes intermolecular distances within 
these complexes and weakens intercluster bonds to so extent that Os clusters can be approximately considered as 
(magnetically) noninteracting units. 

From general point of view, magnetic collapse (disappearance of magnetic properties) observed in e-phase may 
result from coupling of 2, 3, or any other number of O2 molecules in a singlet spin state. The tendency of the O2 
molecules to form such multimolecular clusters (consisting of 2, 3, 4 units) was ascertained long ago in the optical 
spectra of a-phase, where two-, three- (at higher temperature) and four-molecule dipole transition bands were directly 
observed and identified 15 . Why, then, 4O2 complex is more favorable^ than, say, dimer, 2O2, or trimer, 3O2? 

One of the possible reason for such a behavior is a weakness of van der Waals intermolecular forces in comparison 
with exchange interaction. As long as exchange interactions are not taken into account, O2 molecules in a-, f3-, 5- 
and e-phases can be considered as noninteracting solid spheres packed in the most compact way, i.e. in a regular 
triangular lattice, within the basal planei£(see also RefsUIQ). Singlet complexes concatenated by the exchange forces 
and decoupled from each other may also be considered as noninteracting (or weakly interacting) solids. Dimcrs 
themselves^ are highly anisotropic, formation of the decoupled pairs should produce additional distortion of crystal 
lattice (see FigQ]a), so, O4 complexes seem to be unstable with respect to formation of herribone chains. In contrary, 
3O2, 402 and 7O2 complexes may be invariant with respect to rotation around 3-rd, 4-th, or, correspondingly, 6th 
order symmetry axis, and so are isotropic in the basal plane. In turn, a hexagonal plane can be completed by the 
regular triangulars (Fig. QJ>), 60°-angled diamonds (Fig. [J), or hexagons that by appropriate deformations may be 
transformed into highly symmetric 71-O2 units. Lattice distortions shown in Figs. [1] and [2] by arrows could be classified 
(see Tabled according to symmetry of different optical modest. 

It is quite obvious that formation of trimers and sestets should be accompanied by isotropic contraction of inter- 
atomic distances within the complex, while formation of quaters (see below) is related with anisotropic (shear) 
deformation of corresponding rhombus. The second process seems to be energetically more favorable because shear 
modes are usually much more soft compared to isomorphic striction. 

So, formation of quaters may be induced by increase of AFM exchange coupling at high pressure and softening of 
appropriate optical mode. 

II. MODEL 
A. Order parameter 

Phenomcnological description of e-02 as a phase in a series of/3— > a —t 5 — > e transition is not so straightforward 
as of the other magnetic phases. According to phase diagram^, e-phase may be obtained from both 8- and /3-phases 
whose Brave lattices belong to the different space groups, so, what phase should be considered as a parent phase? 
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a) b) 



Figure 1: Covering of hexagonal plane with a) pairs and b) triples of O2 molecules. Arrows show directions of the molecular 
shift in a corresponding optical mode, red parallelogram is a unit cell of the superstructure. 




Figure 2: Covering of hexagonal plane with quarters (outlined with ellipses) of O2 molecules. Arrows show directions of the 
molecular shift in a corresponding optical mode, red rectangle is a unit cell of the superstructure. 

Three-dimensional space group of e-02 coincides with that of a-C>2, though in the phase diagram both phases are 
separated with the high symmetry ^-phase. What is more, a-, S- and e-phases are described by the same symmetry 
group within the basal plane. This fact makes questionable the choice of the components of deformation tensor as an 
order parameter of fe-transition. 

The easiest way to overcome these difficulties is to accept that all the magnetic phases, including £-02, originate 
from a virtual nonmagnetic phase viewed as a stack of regular triangular planes. This assumption is based on the 
following facts. 
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Table I: Wave vectors and polarization of optical modes coupled with different singlet states. Representatives of stars k are 
given according to Kovalev's notations^. Lattice vectors ai, a2 and reciprocal lattice vectors bi, b2 are attributed to the 
hexagonal pra-phase. 



Number of O2 mol. 


Wave vector 


Polarization vector 


2-dimer 


kia = bi/2 


&i + a 2 


3-trimer 


k u = (bi+b 2 )/3 


y r 3(a 1 + a 2 ) + i(ai - a 2 ) 


4-quater 


k 7 = br/2 + ba/4, 
bi/2-b a /4 


u x y/3(ai + a 2 ) + iu y (a 1 - a 2 ) 
iu x V3(a.i + a 2 ) — u v (ai — a 2 ) 



i) Crystal lattices of (3-, a- and <5-phases can be thought of as the different modifications of the same hexagonal 
(space group 6/mmm) pra-phase in which the neighboring close-packed planes are shifted in [1100] direction 10 . 

ii) Though the O2 lattice in e-phase is strongly distorted within the basal plane, as compared to lattice of a-phase, 
an angle between the bonds connecting the molecules in neighboring Os clusters remains approximately equal 
to 60° in a wide interval of pressures, as seen from the experiment^. 

From this point of view, the structural order parameter of pra-phase - e-phasc transition can be represented by the 
amplitudes u x , u y of the optical mode 



u(n) 



ibin/2 



u x cos ■ 



b 2 n 



b 2 n 
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(1) 



where vector n denotes position of a molecule within the basal plane 2 -. Macroscopic description of the magnetic state 
of e-phase may be done with the use of spin-spin correlation functions. Discussion of this question is beyond the scope 
of the present paper. 

It is interesting to note that the structural order parameter in the sequence of f3 — > a — > <5-phase transitions (i.e. 
a function of mutual shift of the neighboring close-packed planes in [1100] direction calculated with respect to the 
initial non-shifted hexagonal stacking) is symmetry related to the transverse acoustic modes propagating in [0001] 
and [1000] directions (wave vectors parallel to b3 and bi, correspondingly). 



B. Free energy and spin hamiltonian 



Different phases of solid oxygen and inter-phase transitions are described on the basis of phenomenological expression 
for free energy of the crystal. Substantial simplification of the model may be achieved by neglection of interplane 
interactions. This assumption is justified by noticeable difference in the variation of in-plane and interplane distances 
in the course of pressure- induced phase transitions. 

As it was mentioned above, the magnetic and crystal structure of a- and (S-phases is indistinguishable within the 
a6-plane, so, herewith we consider the series of (3 — > a — ► e-transition. Gibbs free energy $ of the crystal is modeled 
as a function of (two-dimensional) phonon amplitude u(k), strain tensor components Ujk, invariant with respect to 
the symmetry group of hexagonal pra-phase, plus magnetic contribution into internal energy E mag : 

^ = \ X ( k ^|u(k,-)| 2 + 3 (U XX + Uyyf + ~[{U XX ~ Uyyf + A^y} + P{U XX + Uyy) 

3 

+ J2 Ag° ) (k J -)|ll(k J )| ! W + Uyy) + A( a h n) (k 7 )[^ (k7 ) _ U 2 y (k 7 )](U XX - Uyy) + £ mag . (2) 

3 

Vectors kj in the above expression denote different wave vectors, classified according to irreducible representations 
of 6/mmm space group, phenomenological constants K(kj) are proportional to the corresponding phonon frequen- 
cies, coefficients A p h(k?) originate from the crystal anharmonicity and describe nontrivial coupling between phonon 
amplitude and crystal lattice parameters. Last term in ^ accounts for the (external) hydrostatic pressure P. 

Magnetic contribution E mas is calculated as an average of spin- hamiltonian TL over a ground state |$) of the crystal: 
£ma g = where 

H = J(rnm)S„S m , (3) 
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and summation is accomplished over the nearest and next to the nearest neighbors separated by distance |r nm |. 
Magnetoelastic part of the internal energy is derived from the expression ([3]) with due account of space dependence 
of the exchange integral J(r nm ). 

Once the ground state of magnetic subsystem is calculated, structure and stability conditions of a phase can be 
determined by minimization of free energy ([2]) with respect to phonon amplitudes and deformation tensor components. 

As it was already mentioned, [3- and a-phases posses a kind of the Neel spin ordering, that can be described 
macroscopically by assigning an average value (S n ) to the spin vector at each site. 




Figure 3: Primitive cell of e-phase. Magnetically coupled quaters of molecules are outlined by quadrates. Red and blue arrows 
show the shift of O2 molecules in the course of phase transition. Vectors ai,2 are the Brave lattice vectors of the hexagonal 
pra-phase. 

A ground state l*^) of e-phase is a true eigen function of spin-hamiltonian ([3]) and is calculated within an assumption 
of magnetically decoupled Og clusters. In other words, |\l/ e ) may be represented as unentangled combination of quaters 
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wave functions \ip a ), that satisfies equation 

#(tat»)|tf e } = E ( intra )|* e ), |* e ) EE [] |Vn), (4) 

n 

with the hamiltonian of intra-cluster interactions written as follows 

11S211 + S2 n S3n + S3 n S4 n + S4 n Si n ] 

n 

+ J(r 13 )S ln S 3n + J(r 24 )S 2n S 4n }. (5) 

Here the vectors n define the positions of Os cluster in a superstructure with lattice vectors &[ = 2&i and a 2 = 2a 2 . 
For the sake of simplicity we use the rectangular "unit cell" which contains two clusters (see FigJ3]). Choice of the unit 
cell corresponds to one of three different domains of a-phase. The positions of the individual 2 molecules (labelled 
with number 1, 2, 3, 4) within the cluster are defined with the basis vectors ±ti ee ±a 4 and ±r 2 = ±a 2 . 

Interaction between the clusters with account of the next-to nearest neighbors is described by operator 7^( mtor ) (see 
Fig|3]for notations) 

-^(intci) _ y |j(r' 12 )[Sl„(S 2 (n-n) + S 4 (n+r 2 )) + S 3n (S 2(n _ T2 ) + S 4 ( n+ri )) 
n 

+S 2n (Si( n+Tl ) + S 3 ( n+T2 )) + S 4n (S 1 ( n _ T2 ) + S 3 ( n _ ri ))] 

+ '^( r 'l3)(SlnS 3 ( n+7 - 1 _ T2 ) + S 3n S 1 ( n _ Tl+T2 )) 

+ J( r 24:) Si n (S 3 ( n+ri+T2 ) + S 3 ( n _ Tl _ T2 )) + S 3n (S!( n+ri+T2 ) + S!( n _ Tl _ T2 )) 

+^( r i3) S 2n (S 4 ( n+Tl _ T2 ) + S 4 ( n _ Tl+T2 )) + S 4n (S 2 ( n+Tl _ T2 ) + S 2 ( n _ Tl+T2 )) 

+ ^( r 24)[S2nS 4 ( n+Tl+T2 ) + S 4n S 2 ( n _ Tl _ r2 ))]| , (6) 

so that spin-hamiltonian ([3]) is represented as a sum: 

^ =? ^(i„tra) +? ^(intor)^ ^ 

In the e-phase the first term in ([7]) is responsible for formation of the ground state, while the second one describes 
contribution that arises from excitations. In a-, S and /3-phases both terms contribute equally into magnetic energy 
of the crystal. 

III. e-PHASE 
A. Magnetic structure 

It was already mentioned that according to experimental data, the units Os form a common singlet state, while 
each 2 molecule possesses spin Sj = 1, j = 1 — 4. So, it is convenient to express spin-state tp of Os cluster in terms 
of the eigen functions |0), | ± 1) of spin operators Sf , where Z is a quantization axis. 

According to general theorem of quantum mechanics, singlet state in such a system has 3 representations (among 
81 basic vectors) with the spin wave functions that could be easily found from the conditions 




singlet 



= 0. (8) 



Obviously, V'singict is also an eigen function of hamiltonian 7Y( intra ). 

Additional simplification of the problem may be achieved by account of permutation symmetry group. All three 
singlet states should have different symmetry with respect to permutations of molecules within the cluster and hence, 
correspond to the different eigen values of operators (see Tabic [TT| 

A = (Si, s 2 ) + (s 2 , s 3 ) + (s 3 , s 4 ) + (s 4 , so 
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Table II: Eigen values of operators Pi, P2, P3 in a singlet subspace, cluster energy (7^( mtra )} /N per molecule for arbitrary 
intermolecular spacing, equilibrium angle ip between intermolecular bonds within the cluster and corresponding equilibrium 
energy E £ per molecule. 



Function 


A 


P 2 ,3 


<?i: (intra) }/iv 




E e 


.(singlet) 


-6 


1 


[J(n 3 ) + J(r 24 ) - 6J(n 2 )]/4 


tt/4 


[J(VSo) - 3J(a)]/2 


, (singlet) 
Y-cxl 


-4 





-J(r 12 ) 


arb. 


-J(a) 


.(singlet) 
W2 





-2 


-[J(r 13 ) + J(r 24 )]/2 


tt/4 


-J(V2a) 



P 2 = (S 1; S 3 ), P 3 = (S 2 ,S 4 ) 
Finally, the singlet wave functions may be written in the following form: 

[2|0000) + |010T) + |1010) 



^ g i e t) = _L ( | 1IlT) + |TiT 1)) + _l 



|0101) + 1 1010) - -(10011) + |0110) + 1 1100) + |1001) + 10011) + |0110) 



11100) + 11001)) + -(|1111) 



11111) 



hud 



11111)) 



(9) 



(10) 



, (singlet) 
Vexl 



1 



2\/3 



|1100) + 10011) + |1100) + 10011) 



hud 



hud 



|1111) - |1111) - |1001) - |0110) - 1 1001) - |0110)] 

1 



.(singlet) 

W2 



|1111) + 1 1111) + |1111) + 11111) + 10000) 



|010T) - 10101) - |1010) - |1010)] 



(11) 



(12) 



In order to find out what of three functions (fT0|) -(fT2 |) describes the ground state of hamiltonian ([7|), we compare 
corresponding eigen values (see Tabldrfl the 4th column). Taking into account antiferromagnetic character of the 
exchange interaction (J(r) > 0), the fact, that J(r) monotonically decreases with intermolecular distance r, and 
geometrical relation r 12 < r 13 < r 24 , one can easily verify that 



(H 



(intra 



gr < (H 



(intra 



/cxl 



< (H 



(intra) 1 



cx2 i 



(13) 



and the required ground state is tpg™ &let ^ (see eq.lfTOj)). We have also implicitly taken into account an obvious fact 
that an average value of (7v mtcr )) in any singlet state is exactly zero. 

It is interesting to compare ^ T ln&lct ^ with the AFM Neel state observed in a-phase, where all the nearest neighbors 
are coupled antiferromagnetically. In terms of 2 spin states it means, that the most preferable combinations are 



11111) and 11111) 



States VS'e? 

. (singlet) 



are orthogonal to a subspace spanned over the "Neel-state" vectors | llll) and 

1 1111). In contrary, ip^- 111 ^ 1 '' belongs to this subspace with 0.4 probability. 

We may also compare eigen values of hamiltonian (0) in a singlet state for different clusters: hypothetical dimer, 
trimer and already described quater. It is obvious that for 20 2 and 30 2 complexes the singlet state is kept by the 
nearest neighbor interactions only, corresponding eigen values of hamiltonian (that could be found without explicit 
expression for wave function) are 



(?i (intra) ) d 

^(int ' 



trim 



= -NJ(r 12 ), 

= -N 



for dimer; 



H J(r 12 ) + ij(ri 3 ) 



for trimer. 



(14) 



where N is the number of 2 molecules. In the case of monotonically decreasing AFM exchange and fixed r 12 < 
fi3 < r 2 4 values 



^ (intra) ^(singlet) < (7i^ intra ^) ( ii r 



^(intra^ (singlet) ^ ^ 



(intr 



< (H 



(intra) \ (singlet) 
/ex2 



(15) 



So, the magnetic energy of the crystal in the state with S — takes on its minimum value when the number of 2 
molecules in a singlet group is at least n =4. 
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B. Distortion of crystal lattice 

It was already mentioned that magnetic interactions in solid oxygen are so strong that they cause large distortion 
of crystal lattice. This effect was observed, e.g., in the course of a/3-transitions where in-plane lattice deformation 
achieved nearly 5%2£. In the e-phase the effect of magnetoelastic interactions is even more pronounced, though very 
unusual, because in this case lattice distortion is produced by magnetic collapse, not by magnetic ordering. 

In-plane structure of e-phase can be considered as a result of two-step distortion^ of the ideal hexagonal basal 
plane of pra-phase (/?-) with the lattice constant ah (see Fig|4j: i) homogeneous deformation which changes scales in 
X and Y direction: 

r i3 = a h(l + u xx ), r$ = V3a h (l + u yy ); (16) 

and ii) inhomogeneous distortion of rhombus formed by the in-cluster molecules 1, 2, 3, 4. In fact, this means that 
the virtual intermediate state (after step i)) has an a-type lattice. 




b) c) 

Figure 4: Two-step distortion of the crystal lattice of hexagonal pra-phase (a): i) homogeneous deformation (b); ii) inhomo- 
geneous distortion of rhombus (c). 

Symmetry condition that the distances between the pairs of molecules jk — 12, 23, 34, and 41 within the cluster 
are equal, rjk = a, makes it possible to introduce very convenient and obvious parametrization using an angle tp 
between the directions to nearest neighbors: ri3 = 2acosip, r2A = 2asinc ( 9, and 

u x = a cos <p — 2 T i3 = a(cos <p — cos ipo), u y — ~^ r 2A ~ a sm f = a ( sm <Po — sm <p): (17) 

where rj^ = 2acosipo, r$ — 2a sin tp$ are interatomic distances in the reference frame. 

The positions of O2 molecules in e-phase are then calculated by minimization of Gibbs' free energy ^ with respect 
to the components of deformation tensor Ujk and angle tp. Last term, -Emag, is the magnetic energy in the singlet 
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ground state, 

N 



£ma g = — [J(r X3 ) + J(r 24 ) - 6 J(r 12 )] . (18) 



Thus, expression ((21) can be rewritten as 

^ = ij(Cll + C 12 ){U XX + Uyy) 2 + ljC'[(U XX ~ Uyy) 2 + AU^y] + P(U XX + Uyy) 



2a sin 



-K"(k 7 ) + X^ \k 7 )(u xx + u y y) - A^ n) (k 7 ) cos(<p + tp ){u xx - u yy) 



N 

+ — [J {2a sin ip) + J(2a cos <p) - 6 J(a)] . (19) 

Minimization conditions {d&/d£j = 0, |S 2 <&/e?£j9£k| > 0, £j = (p,u xx ,u yy ) give rise to the following equation for 

dJ(ri 3 ) dJ(r 24 ) . . 4fif(k 7 )a 

— ^ cosyH sm<p-\ — smup — tpo) = (20) 

ar ar Jy 

Assuming the softening of the optical mode k 7 in the vicinity of phase transition point, so that K (k 7 )a 2 <C 
a\dJ(r)/dr\N , we can neglect last term in (|20)l . Then, equation (j20f has an obvious solution i^ cq = 7r/4 (and 
automatically, ri3 = r^Aj- This means that the four molecules in the ground singlet state are situated in the corners of 
quadrate and it is the exchange interaction within the cluster that keeps the molecules in that state. Such a symmetric 
arrangement of molecules seems to be quite natural in the case when the exchange forces are the strongest interactions 
in the system. Really, in the ground state (I10|) the molecules in the neighboring corners (12, 23, 34, and 41) with high 
probability have opposite spins and thus are attracted to each other, due to antiferromagnetic character of exchange 
forces. In the contrary, the molecules in the opposite corners (pairs 13 and 24) have parallel spins and are therefore 
repulsed. The energy of repulsion is minimized when the average distant between corresponding molecules is as much 
as possible, this can be achieved in a symmetrical combination like quadrate. Small deflection (e.g., 96° and 84° at 
17.6 GPa) from the right angle observed in the experiment^ may be calculated from (f2T)]) with account of contribution 
from the optical mode: 

7T 4V2K(k 7 )asin(7r/4-v3 ) , 01 v 

v = ~a ' (21) 

According to Reffl3l. ipo = arctan(r 2 4'* ) = 53°. Using the most elaborated phenomenological for m 21 i 22 of space 
dependence for 

J{r) = J exp[-a(r-r )+/3(r-r ) 2 ], 2.6 < r < 4.2, (22) 

with J = 60 K, a = 3.5 A^ 1 , (3 = 1.2 A~ 2 , r = 3.1854 A, and taking a = 2.18Ai^^, we get an upper limit 
for Jf(k 7 )/iV < 9.2 K/A 2 , while estimated value of J'{r\3)/a > 80 K/A 2 . Considering if as a stiffness constant of 
intermolecular bonds, one obtains the characteristic frequency 16.2 cm -1 which is much smaller than the frequencies 
of optical modes (> 300 cm -1 ) calculated in Refll4l and the frequency 1,38 cm -1 of Raman mode corresponding to 
the antisymmetric stretching motion of the four 2 molecules, coupled in diagonal pairs 12 . 
Stability condition of the "quadrate" solution 

d 2 J(r) 2K(k 7 ) , , 

+ cos( 7 r/4 - po) > (23) 

is obviously satisfied, because according to Ref@ J{r) is a monotonically decreasing concave function (see e.g. (|22]) ) 
of intermolecular distance, and K (k 7 ) > (from the condition of crystal lattice stability). 

So, even in magnetically neutral state the exchange interactions play a role of a motive force that changes crucially 
an angle (p between intermolecular bonds. 

Analysis of the expression Q19p makes it possible to calculate shear deformation u xx — u yy and isotropic striction 
u xx + u yy within the plane: 



2a 2 A( a h n) (k 7 ) 



c 



cos(tt/4 + cp ), 



^ -[P + 2A(r ) (k 7 ) a 2 sin 2 ^i-^] (24) 



Cll + C12 
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Space dependence of the exchange constant J(r) does not contribute into macroscopic deformation, because Os 
clusters are supposed to be decoupled from each other. So, shear deformation of e-phase is due solely to anharmonicity 

(coupling constant X^\k.r)) of crystal lattice. Isotropic striction u xx + u yy describes relative change of the in-plane 
square. From the pressure dependence of lattice parameters^ we can estimate the in-plane compressive modulus 
en + C12 = 88 GPa. From the value of jump of isotropic striction in the fe-transition point, A(u xx + u yy ) = 0.019 we 

estimate isotropic anharmonicity constant (k 7 ) = 1.1 • 10 5 K/A 2 . 



IV. COMPARISON WITH a- AND 5-PHASES 



In the previous section is was shown that once the singlet ground state is formed, crystal lattice should be distorted 
in a described manner, due to strong exchange interactions and reduced optical phonon frequency. But what about 
the inverse mechanism, can the crystal lattice of J-(a-phase) be unstable with respect to u(k 7 ) distortions? 

To answer this question, we minimize Gibb's potential ^ assuming the presence of the collinear long-range AFM 
order. In this case (S%) = 2, (S„) = ±1 at a site R n , and spin polarization alters from 1 to -1 when shifted through 
the vectors ai, a 2 . Taking into account the locations of molecules 1-4 (see Figl3j), it can be easily seen that 



<S!S 2 ) = (S 2 S 3 ) = <S 3 S 4 ) = (S 4 Si) = -1, (S1S3) - (s 2 s 4 ) = 1 
Substituting these values into ([7]) we obtain expression for the magnetic energy Eafm of AFM state: 

N 



(25) 



E, 



may 



E 



AFM 



77iclustcr 1 7^int 
^ AFM + ^ AFM 



[J(r 13 ) + J(r 24 ) - 4J(r 12 )] 



+ j [J(r[ 3 ) + J(r 24 ) + 2(J(r'/ 3 ) + J(r 2 ' 4 )) - AJ(r' 12 )\ 



(26) 



where the first term describes interactions inside the cluster and the second one is responsible for interaction energy. 

Using the same parametrization of shift components u x , u y as in (|1T|) we express all the intermolecular distances 
(see FigE]) in the expression (f2"6"| in terms of ip, a as follows: 



r' 12 = a\J (2 cosyo — cos p) 2 + (2 sin^o — simp) 2 , 



(27) 



r' 13 = 2a(2 cos po — cos tp), r" 3 = 2a\/ cos 2 (po + (sin ipo — sin ip) 2 , 
2a(2sin(,!3o — sirn^), r 2 ' 4 = 2aJ sin 2 ipo + (cosi^o — cosy) 2 . 



' 24 



Analysis of the expressions ([2]), ([26]) shows that the conditions of minimum d$/d(p = 0, d 2 $>/dtp 2 > for AFM state 
are satisfied for tp = cpo, r' 12 = ri 2 , r' 13 



' _ r " — T .(°) r ' 
'13 — '13 ' '24 



r 24 — r 24 , as can be seen from the following relations 



JT = Na 
o<p 



J'(r[ 2 ) 



+ J'(r 13 ) 



2 sm(p a - ip) 
\/5 - Acos(p - ip ) 

cos p(sin ip — sin <po) 
a/cos 2 ipo + (sin^o — simp) 2 



J'(r M ) 



sin y(cos ipo — cos ip) 
\J sin 2 ipo + (cos ipo — cos ip) 2 



(28) 



2a 2 A'(k 7 ) sin(y - <p ) = 



— j — Na [J'(r 24 ) sin tpo + J'(r 13 ) cos^o — 2J'(r' 12 )] + 2a 2 A'(k 7 ) > 0. (29) 
dip* 

Inequality (|29| is obviously satisfied due to the already mentioned fact that the exchange integral is a positive and 
monotonically decreasing function of intermolecular distance. 

So, in AFM state the crystal lattice is stable with respect to distortions even in the case of vanishingly small stiffness 
A"(k 7 ). As in the case of e-phase, it is the exchange forces that keep the lattice from distortion. In the phase with 
the long-range magnetic ordering the values of the "exchange bonds" pulling the 2 molecules in opposite directions 
are equal (compare with e-phase, FigH]), and this impedes nonsymmetrical distortion of Os-rhombuses. 

It is also instructive to compare the magnetic energies of AFM 



Eafm = N[J(2acosipo) + J(2asimpo) — 2.7(a)], 



(30) 
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and e-phases 

N r- 

E e = —[J(V2a)-3J(a)}. (31) 

It is obvious that in a nondeformed hexagonal lattice (<po — 60°) E s < Eafm for any value of a. This means that 
the AFM state of a and <5-phases is stabilized by the long-range elastic forces that produce homogeneous deformation 
(striction) of crystal lattice, as it was shown in Refill 



V. CONCLUSIONS 



In summary, we have calculated the wave functions of singlet state implemented on the 4(02) cluster and found 
that the exchange energy of the ground state (|10j) is lower than that in another singlet states implemented on dimers, 
trimers and quaters. Interactions between next-to-nearest neighbors (i.e., between O2 molecules located in the opposite 
corners of rhombuses) plays an important role in stabilization of the magnetic and crystal structure in the ground 
state. 

We have shown that the observed distortion of crystal lattice and formation of 4(02) quadrates in e-phase can 
be explained by strong magnetoelastic contribution into exchange energy along with the softening of u(kj) optical 
phonon mode. Stability of the distorted lattice in e-phase is then due to antiferromagnetic character of exchange 
interaction in solid oxygen. 

The same magnetoelastic forces ensure stability of the AFM long-range phases (a, 5) with respect to inhomogeneous 
distortion of crystal lattice even in the case when stiffness constant of u(k7)-mode is vanishingly small. 
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In the absence of magnetic interactions a5-transition can be classified as a Ilnd order, according to Liftshits' criterium. 
Disappearence of the magnet ic p roperties of solid oxygen was experimentally proved in Reffl9l. 

In the very recent paper Refl23l the authors claim that according to their quantum chemical calculations a "rhomboid Os 
structure of D2/1 symmetry is a stable species, that is, a local energy minimum" and point out that the ground state of a 
rhomboid should be a singlet spin state. On the basis of density functional calculations the authors also predict stability of 
the trimer structure S6 as corresponding to a local energy minimum on the potential hypersurface. 

In contrast to homogeneous deformations, an optical mode removes degeneracy of the nearest neighbors intermolecular 
distance. 

Order parameter of pra-to-/3-, a- and 5- phase transitions is coupled with a transverse acoustic mode in [0001] direction. 
In approximation of linear elasticity the sequence of steps is immaterial. An accepted sequence is convenient from method- 
ological point of view. 
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